    local controls lnsize age1900 fxassratio_MEAN  
    local depvar DR 
    local varlist1 denom  
    local varlist2 Dstocksize_top3rd  `controls'
    local varlist3 Dstocksize_low3rd  `controls'
    local varlist4 Dstocksize_low3rd  `controls' founderfirm
    local varlist5 Dstocksize_low3rd  lnsize leverage_MEAN
    
    logit `depvar' `varlist1' , nolog robust     
    margins, at(denom=(237.801 2362.199))  asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  

    logit `depvar' `varlist2' , nolog robust     
    margins, dydx(`varlist2')  asobserved   

    qui logit `depvar' `varlist2' , robust     
    margins, at(lnsize=(.754 1.377)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  

    qui logit `depvar' `varlist2' , robust     
    margins, at(age1900=(37.43 64.26)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  

    qui logit `depvar' `varlist2' ,  robust     
    margins, at(fxassratio_MEAN=(.312 .488))  asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  
 
 
    logit `depvar' `varlist3' , nolog robust     
    margins, dydx(`varlist3')  asobserved   

    qui logit `depvar' `varlist3' , robust     
    margins, at(lnsize=(.754 1.377)) asobserved post             /* one std dev around the mean */
     di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  

    qui logit `depvar' `varlist3' , robust     
    margins, at(age1900=(37.43 64.26)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  

    qui logit `depvar' `varlist3' ,  robust     
    margins, at(fxassratio_MEAN=(.312 .488))  asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  

    logit `depvar' `varlist4' , nolog robust     
    margins, dydx(`varlist4')  asobserved   

    qui logit `depvar' `varlist4' , robust     
    margins, at(lnsize=(.754 1.377)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  

    qui logit `depvar' `varlist4' , robust     
    margins, at(age1900=(37.43 64.26)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  

    qui logit `depvar' `varlist4' ,  robust     
    margins, at(fxassratio_MEAN=(.312 .488))  asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  

    logit `depvar' `varlist5' , nolog robust     
    margins, dydx(`varlist5')  asobserved   

    qui logit `depvar' `varlist5' , robust     
    margins, at(lnsize=(.754 1.377)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  

    qui ologit `depvar' `varlist5' , robust     
    margins, at(leverage_MEAN=(.297 .491)) predict(pr outcome(1)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]  
    test _b[1._at] = _b[2._at]  
